Anisotropic higher derivative gravity and inflationary universe 
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Stability analysis of the Kantowski-Sachs type universe in pure higher derivative gravity theory is 
studied in details. The non-redundant generalized Friedmann equation of the system is derived by 
introducing a reduced one dimensional generalized KS type action. This method greatly reduces the 
labor in deriving field equations of any complicate models. Existence and stability of inflationary 
solution in the presence of higher derivative terms are also studied in details. Implications to the 
choice of physical theories are discussed in details in this paper. 

PACS numbers: 98.80.-k, 04.50.+h 



I. INTRODUCTION 



Inflationary theory is a nice resolution for the the flatness, monopole, and horizon problems of our present universe 
described by the standard big bang cosmology P|. Our universe is homogeneous and isotropic to a very high degree 
of precision 0, Q . Such a universe can be described by the well known Friedmann- Robertson- Walker (FRW) metric 
i. 

One expects that gravitational physics could be different from the standard Einstein models near the Planck scale 
CN ' [E 111' For example, quantum gravity or string corrections could lead to interesting cosmological applicastions0- 
Indeed, some investigations have already addressed the possibility of deriving inflation from higher order gravitational 
corrections HHSEl- 

In particular, a general analysis of the stability condition for a variety of pure higher derivative gravity theories 
could be very useful in many respects. In fact, it has been shown that a stability condition should hold for any 
■ potential candidate of inflationary universe in the flat Friedmann- Robertson- Walker (FRW) space ^ simple 

method has been shown useful in such analysis. 

In addition, there is no particular reason why our universe is initially isotropic to such a high degree of precision. 
. Even anisotropy can be smoothed out by the proposed inflationary process, it is also interesting to study the stability 
' of the FRW space during the post-inflationary epoch. Nonetheless, it is interesting to study the cases where our 
^ ' universe starts out from an initially anisotropic universe. As a result, our universe is expected to evolve from certain 
I anisotropic universe to an stable isotropic universe. Indeed, it has been shown that there exists such kind of anisotropic 
O i' solution for a NS-NS model with a metric field, a dilaton and an axion field This inflationary solution is also 
D , shown to be stable against small field perturbations 1^. Note also that stability analysis has also been studied in 

' various fields of interest [lA J4j . 
J> ', Higher derivative terms should be important for physics near the Planck scale [lOlIT^. For example, higher order 
corrections from quantum gravity or string theory have been considered in the studies during the inflationary phase 
[l5l |. Higher derivative terms also arise as the quantum corrections to the matter fields Moreover, the stability 
analysis of the pure higher derivative gravity models has been shown in Ref. (l0l| . It is then interesting to study the 
implication of this stability analysis in different models. 

Recently, there are growing interests in the study of Kantowski-Sachs (KS) type anisotropic universes [T^ IT^ ITsjl . 
Hence we will try to study the existence and stability of an inflationary de Sitter final state in the presence of higher 
derivative gravity theory on the Kantowski-Sachs spaces. In particular, it will be applied to study a large class of 
pure gravity models with inflationary KS/FRW solutions in this paper. Any KS type solution that leads itself to an 
asymptotic FRW metric at time infinity will be referred to as the KS/FRW solution in this paper for convenience. 

It will be shown that the existence of a stable de Sitter background is closely related to the choices of the coupling 
constants. We will try to generalize the work in Ref. 0, |23| in order to obtain a general and model-independent 
formula for the non-redundant field equations in the Kantowski-Sachs (KS) type anisotropic space. This equation can 
be applied to provide an alternative and simplified method to obtain the stability conditions in pure gravity theories. 
In fact, this general and model- independent formula for the non-redundant field equations is also very useful in many 
area of interests. 

We will first derive a stability equation which turns out to be identical to the stability equation for the existence 
of the inflationary de Sitter solution discussed in Ref. '20]. Note that an inflationary de Sitter solution in pure 
gravity models is expected to have one stable mode and one unstable mode for the system to undergo inflation with 
the help of the stable mode. Consequently, the inflationary era will come to an end once the unstable mode takes 
over after a brief period of inflationary expansion. The method developed in Ref. [Toll20| was shown to be a helpful 
way in choosing physically acceptable model for our universe. Our result indicates, however, that the unstable mode 
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will also tamper the stability of the isotropic space. To be more specific, if the model has an unstable mode for the de 
Sitter background perturbation with respect to isotropic perturbation, this unstable mode will also be unstable with 
respect to any anisotropic perturbations. 

In particular, we will show in this paper that the roles played by the higher derivative terms are dramatically 
different in the inflationary phase of our physical universe. First of all, third order term will be shown to determine 
the expansion rate Hq for the inflationary de Sitter space. The quadratic terms will be shown to have nothing to do 
with the expansion rate of the background de Sitter space. They will however affect the stability condition of the de 
Sitter phase. Their roles played in the existence and stability condition of the evolution of the de Sitter space are 
dramatically different. 



II. NON-REDUNDANT FIELD EQUATION AND BIANCHI IDENTITY IN KS SPACE 

Given the metric of the following form: 

ds^ = -dt^ + c^{t)dr^ + a^{t){d^e + f{e)dif^) (1) 

with f{d) = sinh 0, sin^) denoting the flat, open and close anisotropic space known as Kantowski-Sachs type 
anisotropic spaces. To be more specific, Bianchi I (BI), III(BIII), and Kantowski-Sachs (KS) space corresponds to 
the flat, open and closed model respectively. One can instead write the metric as 

ds" = ~de + a\t){^^ + r^de^) + al{t)dz^ (2) 

with r, 0, and z read as the polar coordinates and z coordinate for convenience. One writes it this way in order 
to make the comparison with the FRW metric easier. Note that k = 0,1,-1 stands for the flat, open and closed 
universes similar to the FRW space. 

In addition, one can restore the gtt component lP'{t) for the purpose of deriving the non-redundant field equation 
associated with Gtt that will be shown shortly. As a result, one has 

ds2 ^ -b^(t)dt^ + a^(t){ + r^de^) + al{t)dz^ (3) 

1 — fcr^ 

One can show that all non-vanishing spin connections read 
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(5) 
(6) 

(7) 

(8) 

(9) 

Here B = 1/6^ and Hi = {a/a, a/a, dz/az) = {Hi, H2 — Hi, Hz) for r, 6, and z component respectively. One can also 
define = V/^,^ and F'^ = F^J^g''" for convenience. As a result, one can compute all their non-vanishing components, 
with 6=1 being reset to unity: 

F* = -iH (10) 

r'- - A (11) 
Ft = 3H (12) 

= r{l-kr^)- ^^^^ 

Writing iJ^^ = G^^ — T)„y, one can show that D^H'^" = follows from the Bianchi identity D^G^^'^ = and the 
energy momentum conservation D^T^'^ — for any energy momentum tensor T^'^ coupled to the system. With the 
metric 0, one can show that the r component of the equation D^H^'^ — implies that 

i/; = H%. (14) 
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The result says that any matter coupled to the system has the property that T\. = T^g- In addition, the the equations 
Dfj^H^^ — and D^H^^'^ — both vanish identically all by itself irrelevant to the form of the energy momentum 
tensor. More interesting information comes from the t component of this equation. It says: 

{dt + iH)H\ = 2HiH\. + H,H''^. (15) 

This equation implies that (i) H\ = implies that H"^^ = H\ = and (ii) i/"^ — li\ = only implies that 
{dt + 3H)H\ = 0. Case (ii) can be solved to give H\ = constant exp[— a^Oz] which approaches zero when a^Uz oo. 
Therefore, for the anisotropic system one is considering here, the metric contains two independent variables a and 
while the Einstein field equations have three non-vanishing components: H\ — 0, H^^ — H^g ~ and H\ — 0. The 
Bianchi identity implies that the tt component is not redundant which needs to be reserved for complete analysis. 
One can freely ignore one of the rr or zz components without affecting the final result of the system. In short, the 
Hi — equation, known as the generalized Friedman equation, is a non-redundant field equation as compared to the 
= and = equations. 

In practice, the field equations of the system may be very complicate. One of the resolution is to reduce the 
Lagrangian of the system from a functional of the metric gfj,u, C{gfiu), to a simpler function of a{t) and az{t), namely 
L{t) = a?'azL(gfiu{a{t), az(t))) . The resulting equation of motion should be able to be reconstructed from the variation 
of the reduced Lagrangian L(t) with respect to the variable a and a^. The result is, however, incomplete because, the 
variation of a and ttz are related to the variation of Qrr and gzz respectively. One can never derive the field equation 
for gtt without restoring the variable h{t) in advance. This is the reason why one wishes to introduce the metric Q 
such that the reduced Lagrangian Ht) = a'^azC{g^^{b{t),a{t),az{t))) contains the non-redundant information of the 
= equation. One can reset b = 1 after the variation of b{t) has been done. This will reproduce the wanted and 
non-redundant Friedman equation. 

After some algebra, one can also compute all non-vanishing components of the curvature tensor: 

- [^BH,+BiH, + Hf)W, (16) 
R'\i = BH,H,e'^'"eki,n + \e'''ekiz. (17) 



Given a Lagrangian L = ^/gC = L{b{t), a{{t), az{t)) one can show that 



L = ^C{R'\^,R\,) = ^C{H,,H,,a') (18) 



The variational equations for can be shown to be: 

£ + Hi{-^ + 3H)V = H,L,+HiV (19) 

C + i^+iHfL' = {^+3H)L,-a^^il-d^z) (20) 
at at oa^ 

Here Li = S£/5Hi, U = 5C/5Hi, and iH = Hi. For simplicity, we will write C as, L from now on in this paper. 
As a result, the field equations can be written in a more comprehensive form: 

L + H,{j^+3H)V = H,L,+H,V (21) 

L+{^ + 3HfL' - {^ + 3H)L,~a^^{l-S,z) (22) 
dt dt oa'^ 



III. FRW SPACE AS A STABLE FINAL STATE 



For simplicity, one will start with Einstein-Hilbert (EH) action and study its evolutionary process from an anisotropic 
space. It is known that our final universe is isotropic to a very high precision. Therefore, any physical model should 
carry any physical universe from initial anisotropic space to an isotropic space as a final destination. Since lowest 
order Einstein-Hilbert action is the most well-known popular model, one expects such isotropilized process should 
be realized in this model. Any acceptable higher order terms being considered as corrections around this stable EH 
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background should not affect its intention evolving toward the FRW space. Therefore, one will start with a simple 
EH action with a cosmological constant term given by: 



One can show directly that the reduced Lagrangian L is given by 



2k 



2A. 



Therefore, one can show that the Friedniann equation (|21|l and z-equation (|22|l take the following form: 



Hf + 2HiH, + ^ = A, 
2Hi + iHl + 4 = A- 
For convenience, one can also use equations (|25I26|) to derive the following equation: 

Hi + H\ = HiHz- 

Eq. (|27|) can be shown to give 

Qz = kod 

for some integration constant ko. Note also that Eq. (|26|l can also be integrated as 

rr2, k _ A fcl 



(23) 

(24) 

(25) 
(26) 

(27) 
(28) 

(29) 



for some integration constant ki. For the case k = 0, one can show that A = %HqA/A if we write A = c?!"^ . Hence 
one can integrate this equation to obtain 



a{t) = a(0) 



exp[3i?oV2] + fc2 exp[-3i/oV2] 



-1 2/3 



l + fc2 



(30) 



for come constant k-^. Here — A/3 denotes the expansion factor. For the case where fc 7^ 0, a special solution with 
/ci = can be found to be 



a = a\ 



exp[i7oi] 



exp[— i^ot] 



(31) 



In fact, one can show that the evolutionary properties of these solutions can be obtained without knowing the exact 
solutions. Indeed, one will show in a moment that the inflationary solution will try to evolve to an isotropic FRW 
space as i 00. In addition, one can also show that these solution will remain isotropic from an stability analysis. 
One will study the evolution of physical universe from anisotropic initial state to isotropic final state. Therefore, one 
will write above equations in terms of the following variables: 



m = 2Ex + Hz 
A = (i/i - Hz) 



V 



(32) 
(33) 



with T,, V = a^ttz, and A the total expansion rate, 3- volume, and the deviation function respectively. One expects 
A — > as the physical universe evolves toward an isotropic final state. Indeed, one can show that the field equations 
can be written as 



A + 3HA = -- 



(34) 



Eq. (|34(l can be rewritten as 



dt 



(VA) 



-kaz- 



(35) 
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Hence one has 

A = A(0)^_^i^i(^ (36) 
V az 

If az eventually expands proportional to exp[Aoi], one can show that az{t')dt' /a^ — ^ I/Aq, a small constant, as 
t ^ oo. Since the scaling factor a^(0)/a^ for an expanding a solution, this equation implies that A ^ as 
t oo. Hence one shows that the field equations of the EH action will definitely take the anisotropic universe, either 
one of the KS type spaces, to the final FRW universe as t ^ oo. 

One can also show that the final isotropic FRW universe is stable against any small perturbations Hi = Ho + SHi 
and Hz = Hq + SHz- For convenience, one can also use equations H25I26|I to derive the following equation: 

Hz+Hl + 2HiHz = A. (37) 
Applying this perturbations to Eq.s l(77|l and ^7\i . one has 

SHi + Ho5Hi -Ho6Hz = 0, (38) 

6Hz + AHqSHz + 2Ho6Hi = 0. (39) 

Adding equations H38I39|I one can derive 

{SHi + ZHqSHi) + {5Hz + SHqSHz) = 0. (40) 

Hence one has SHi + SHz = constant/a'^ ^ as a — > oo. Hence any physical perturbation against the the FRW 
background would imply SHi —SHz = SH. Hence, Eq. implies 

SH + 2HoSH = 0. (41) 

This equation can be integrated to obtain the result 

constant 

SHi = > 42 

as t oo. Hence one can again show that the final FRW space is stable against this anisotropic perturbation. Note 
that both SHi asymptotically, while their signs in the order of 0{a^^) are opposite to each other such that 

their sum 6 Hi + SHz 



IV. STABILITY OF HIGHER DERIVATIVE INFLATIONARY SOLUTION 



One can then apply the perturbation, Hi = HiQ + SHi^ to the field equation with Hio the background solution to 
the system. This perturbation will enable one to understand whether the background solution is stable or not. In 
particular, one would like to learn whether a KS — > FRW (KS/FRW) type evolutionary solution is stable or not. 

Note that our universe could start out anisotropic even evidences indicate that our universe is isotropic to a very 
high degree of precision in the post inflationary era. Therefore, one expects that any physical model should admit a 
stable KS/FRW solution. In particular, one will be interested in a de Sitter (dS) background solution with HiQ = Hq 
for some constant Hubble expansion parameter. One will denote such solution as KS — *■ de Sitter (KS/dS) type 
inflationary solution. 

One can show that any FRW inflationary solutions with a stable mode and an unstable mode is a negative result 
to our search for a stable inflationary model. In particular, any FRW inflationary solutions with a stable mode and 
an unstable mode will provide a natural way for the inflationary universe to exit the inflationary phase. Such models 
will, however, also be unstable against the anisotropic perturbations. Therefore, such solution will be harmful for 
the system to settle from anisotropic space to FRW space once the graceful exit process is done. One will show in 
this section that the higher derivative gravity theory one considers here could also accommodate two stable modes 
with appropriately chosen coupling constants. In such case, the inflationary de Sitter solution Hq will also be stable 
against anisotropic perturbations. 

First of all, one can show that the first order perturbation equation from the non-redundant field equation l|21|) . 
with Hi ^ Hq + SHi, gives 



< H.L'^SHj > +3H < H.V^SHj > +S < H,t > +3H < {H,L) + V)SHj > 
+ < H,V > S{3H) H,L,jSHj > + < H,SV > 



(43) 
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with all functions of Hi evaluated at some FRW background with Hi = Hq. The notation < AiBi >= J2i=i z ^iHi 
denotes the summation over i = \ and z for repeated indices. Note that we have absorbed the information of z = 2 
into i — 1 since they contributes equally to the field equations in the KS type spaces. In addition, = 5^L/SHi5Hj 
and similarly for L^- and L'^ with upper index * and lower index j denoting variation with respect to Hi and Hj 
respectively for convenience. In addition, perturbing Eq. (|22|l can also be shown to reproduce the Eq. (|43|l in the 
FRW hmit*Q. 

Once we adopt the de Sitter solution with Hq ~ constant, one has 

< H.V^Hj > +3H < H,V^5Hj > +3H < {H,L] + U)5Hj > + < H^V > S{3H) =< H,L,j5Hj > (44) 

where all field variables are understood to be evaluated at the background de Sitter space where Hi = Hq ~ constant 
for all directions. 

If the inflationary de Sitter solution has one one stable mode and one unstable mode for the system, the unstable 
mode is expected to collapse the de Sitter phase. Then the inflationary era will come to an end once the unstable 
mode takes over. It was shown earlier to be a helpful way to select physically acceptable model for our universe. Our 
result shown here indicates, however, that the unstable mode will also tamper the stability of the isotropic space. 
Indeed, if the model has an unstable mode for the de Sitter perturbation, this unstable mode will also be unstable 
against the anisotropic perturbation. 

For example, one can show that the model 

C = -R + aR" + PR^R''^ + ^R^''^,^ R^l^ R^p^^ (45) 

admits an inflationary solution if 7 > 0. Note that the 7 term is the minimal consistent effective low-energy two-loop 
renormalizable Lagrangian for pure gravity theory (2ll | . In addition, the quadratic terms can be shown to be derivable 
from the matter effect of quantum fields. For simplicity, one can write 

A = Hi+Hl, (46) 

B = Hf + ^, (47) 

C = HiH,, (48) 
D = H, + Hi (49) 

since all curvature tensor components and all field equations will be functions of above combinations. This notation 
will be shown to very convenient in tracking the field equations for any complicate models such as the higher derivative 
models we are working on in this paper. Indeed, one can show that the Lagrangian reads 

L = AA + 2B + AC + 2D + 'ia [AA^ + B^ + 4(7^ + + AAB + SAC + AAD + ABC + 2BD + ACD] 

+2P [3A^ + B^ + 3C^ + D^ + 2AB + 2AC + 2AD + 2BC + 2CD] + 87 [2A^ + B^ + 2C^ + D^] . (50) 

This Lagrangian will reduce to the de Sitter models when we set Hi ^ H in the isotropic limit. Hence one can show 
that the generalized Fricdmann equation (|21|l gives 

^0 - 1/47 (51) 

in the isotropic limit for an inflationary solution with ^ 1. 

Note that the quadratic terms does not contribute to the Friedmann equation in the de Sitter limit with a constant 
inflationary phase. This can be shown to be a general property of the quadratic models. One will show later that 
these terms will, however, affect the stability conditions for the inflationary phase. Indeed, quadratic terms will be 
shown to affect the duration of the inflationary de Sitter phase which will be shown in a moment. 

Note that one can also show that 

< H,L'^ > = 2 < >, (52) 

< H,L\ > = 2 < H,Ll >, (53) 

= 2L\ (54) 

< H,La > = 2 < H,L,^ >, (55) 

in the inflationary de Sitter background with Hq = constant. Therefore, the stability equations (|44|l can be greatly 
simpliflcd. For convenience, one will define the operator T) as 

VSH =< H,V^ > SH + 3H < H,V^ > SH + 3H < H,L\ + > 5H + 2 < H,V > SH- < H,La > SH. (56) 
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As a result, one can show that the stabiHty equation (|44|l reads 

V{SHi+SH,) ^0 (57) 

with i = 1 and i — z component summed over. Since, one can perturb the filed Hi from any directions, stability 
conditions must also hold for perturbation from any directions. Hence, one expects the stability conditions must hold 
and be identical from the perturbations in each directions. Indeed, above result shows that both 5Hi and 5Hz follows 
the same condition. Hence one needs to solve the following stability condition VSHj — for both i = 1 and i — z. 
Hence one has 

< HiV^ > 5Hj +3H < HiV^ > 5Hj + 3H < HiL\ + > 5Hj + 2 < HiV > SHj =< H,L,i > SHj (58) 

to see if the system is stable or not against any small perturbation with respect to the de Sitter background. 
In addition, the stability equation H58|) for SHi can be shown to be 

[6a + 2(3+ l2jH^] \SHi + SHaSHi] + (1 - l2-fH^)SHi = (59) 

for such KS/dS solutions. Hence one has 

SH, = exp[^^(l + S,)] + d, exp[^^(l - S,)] (60) 

with 



^1 = y^l + 8/[27 + 9(6a + 2/3)iJ2] (61) 

and some arbitrary constants Ci,di to be determined by the initial perturbations. It is easy to see that any small 
perturbation SHi will be stable against the the de Sitter back ground if both modes characterized by the exponents 

A± = -[3Hot/2][l±Si] (62) 

are all negative. This will happen if 5i < 1. In such case, the inflationary de Sitter space will remain a stable 
background as the universe evolves. 

On the other hand, one would have a stable mode and an unstable mode ii Si > 1. This indicates that this model 
admits one stable mode and one unstable mode following the stability equation (|58|l for the inflationary de Sitter 
solution. It is shown to be a positive sign for an inflationary model that is capable of resolving the graceful exit 
problem in a natural manner. 

Indeed, one expects any unstable mode for a model to be of the form SHi ~ exp[lHot], to the lowest order in HqI, 
in a de Sitter background with I some constant characterizing the stability property of the model. In such models, 
the inflationary phase will only remain stable for a period of the order At ^ I/IHq. The inflationary phase will start 
to collapse after this period of time. This means that the de Sitter background fails to be a good approximation when 
t > At. 

As a result, the anisotropy will also grow according to SHi SHi exp[lHoAt] with SHi denoting the initial 
perturbation. Hence this model will have problem remaining isotropic for a long period of time. Therefore, pure 
gravity model of this sort will not solve the graceful exit problem. One will need, for example, the help of certain 
scalar field to end the inflation in a consistent way. The unstable mode gives us, however, a hope that small anisotropy 
observed today can be generated by the initial inflationary instability for models with appropriate factor I. 

The result shown in this paper shows that the roles played by the higher derivative terms are dramatically different 
in the inflationary phase of our physical universe. First of all, third order term characterized by the coupling constant 
7 will determine the expansion rate Hq, given by Eq. (|51|l . for the inflationary de Sitter space. The quadratic terms 
characterized by the coefficients a and (3 will not affect the expansion rate of the background de Sitter space. They 
will however affect the stability condition of the de Sitter phase depending on the sign of the characteristic function 
A-I-. Both the third order term and quadratic terms are closely related to the quantum corrections of the quantum 
fields [TM I2H. Their roles played in the existence and stability condition of the evolution of the de Sitter space are 
dramatically different. They are however equivalently important in the higher derivative models. 



V. CONCLUSION 



We have tried to generalize the work in Ref. T^, ^2d\ in order to obtain a general and model-independent formula for 
the non-redundant field equations in the Kantowski-Sachs (KS) type anisotropic space. This equation can be applied 
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to provide an alternative and simplified method to obtain the stability conditions in pure gravity theories. In fact, 
this general and model-independent formula for the non-redundant field equations is also very useful in many area of 
interests. 

It is also shown that the existence of a stable de Sitter background is closely related to the choices of the coupling 
constants. We first derive a stability equation which turns out to be identical to the stability equation for the existence 
of the inflationary de Sitter solution discussed in Ref . |0, . 

If the inflationary de Sitter solution in the pure gravity theory has one one stable mode and one unstable mode 
for the system, the unstable mode is expected to collapse the de Sitter phase. Consequently, the inflationary era 
will come to an end once the unstable mode takes over after a brief period of inflationary expansion. The method 
developed in Ref. 0,1131 is useful way in choosing physically acceptable model for our universe. Our result indicates, 
however, that the unstable mode will also tamper the stability of the isotropic space. 

To be more specific, if the model has an unstable mode for the de Sitter background perturbation with respect 
to isotropic perturbation, this unstable mode will also be unstable with respect to any anisotropic perturbations. 
In particular, we have shown in this paper that the roles played by the higher derivative terms are dramatically 
different in the inflationary phase of our physical universe. First of all, third order term is shown to determine the 
expansion rate Hq for the inflationary de Sitter space. The quadratic terms will be shown to have nothing to do 
with the expansion rate of the background de Sitter space. They will however affect the stability condition of the de 
Sitter phase. Their roles played in the existence and stability condition of the evolution of the de Sitter space are 
dramatically different. 

In short, the result of this paper shows that graceful exit and stability of any de Sitter model can not work along 
in a naive way. The physics behind the inflationary de Sitter models appears to be much more complicate than one 
may expect. In another words, the phase transition during and after the inflationary phase deserves more attention 
and requires extraordinary care in order to resolve the problem lying ahead. 

In addition, the result shown in this paper shows that the roles played by the higher derivative terms are dramatically 
different in the inflationary phase of our physical universe. First of all, third order term characterized by the coupling 
constant 7 will determine the expansion rate Hq, given by Eq. H51II . for the inflationary de Sitter space. The quadratic 
terms characterized by the coefficients a and /3 will not affect the expansion rate of the background de Sitter space. 
They will however affect the stability condition of the de Sitter phase depending on the sign of the characteristic 
function A± . Both the third order term and quadratic terms are closely related to the quantum corrections of the 
quantum fields 15, 21]. Their roles played in the existence and stability condition of the evolution of the de Sitter 
space are dramatically different. They are however equivalently important in the higher derivative models. 
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